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ABSTRACT 

We present new observations of the prominence system on the K3 dwarf Speedy Mic (BO 
Mic, HD 197890). Using an improved technique to track the absorption features in Ha we 
find a very active prominence system with approximately ten prominences on the observable 
hemisphere per rotation. From a total of 25 prominences we find an average axial distance of 
(2.85 ± 0.54)_R* which is twice the co-rotation radius above the stellar surface. We discuss 
the consequences of these observations on the nature of the supporting magnetic structures. 
Two consecutive nights, with complete phase coverage, combined with a further night after 
a three night gap allow us to study the evolution of the prominence system on two different 
timescales. Several of the prominences have counterparts at similar phases on consecutive 
nights. During this interval many prominences show evidence for evolution in their heights 
and phases of observation. Five nights (13 rotation cycles) later we recover many prominences 
at approximately the same phases. Whilst individual prominences change axial distances or 
appear/re-appear from night-to-night the underlying prominence supporting structures appear 
to be stable over as many as 13 stellar rotations. 
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1 INTRODUCTION 

Studies of young rapidly rotating cool stars have revealed 
many solar ana logues. Advances in Doppler imaging techniques 
jVoet & PenrodJ983.) have now shown a number of highly spot- 
ted st ars indicating efficient magnetic dynamo activity iStrassmeieil 
l2002h . More recently Zeeman-Doppler Imaging (ZDI) has pro- 
vided the first magnetic surfac e maps of these active stars 
jPonati & Collier Cameronll997^ . In order to fully understand the 
field topology of rapidly rotating stars, we need to examine the re- 
lationship between surface and coronal magnetic fields. Potential 
field extrapolations have made use of magnetic images to model 
the coronal magnetic structure I Jardine et al. 2002). This enables a 
comparison to be made with the hot extended active stellar coronae 
inferred from X-ray observations. Here we present observations of 
stellar prominences which give us a direct probe into the coronal 
magnetic structure and compare our results in the context of field 
configurations. 

Stellar prominences are observed as rapid spectroscopic varia- 
tions in the stellar Balmer lines. Transient absorption features pass- 
ing through the Dopp ler broadened Ha line of rapidly rotatin g stars 
were first reported bv lColIier Cameron & RobinsonHl989ah on the 
rapidly rotating KOV dwarf AB Doradus (AB Dor). These were 
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found to be cool dense condensations of material in enforced co- 
rotation with the stellar surface and so dubbed "prominences". 

From well sampled Ha timeseries one can infer the spatial 
distribution of the prominence structures from the rate at which 
absorption features are seen to drift through the profile. Such ob- 
servations have been carried out by different authors on around 
a dozen stars over the last two decades. For a recent review of 



prominence observations see 'Cameron et aL'^OOS*). AB Dor how- 
ever, remains the most studied object to date. Prominences have 


consistently been reported { 
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range of distances from the stellar rotation axis from 2 to 8 i?*, 
with a concentration at the co-rotation radius (2.7i?* on AB Dor). 

As their solar namesake suggests these stellar prominences are 
confined and supported by magnetic loop structures. This is an im- 
mediate result of finding prominences beyond the co-rotation ra- 
dius where a magnetic tension term is needed to balance the out- 
ward centrifugal force. Therefore observations of prominence loca- 
tions and evolution provide a unique probe to the stellar magnetic 
field out to distances of many stellar radii into the corona. 

Here we present new observations of the prominence system 
on one of the fastest and brightest young rapidly rotating stars in the 
solar neighbourhood. Speedy Mic (BO Mic, HD 197890) is a K3V 
dwarf with an equatorial rotation velocity wsini~ 132kms~^ 
jBameslEool . Recent studies of X-ray variability have shown 
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Figure 1. An example observed Ha profile (solid line) shows an absorption 
transient. Also plotted is the average profile (dotted line) for comparison. 
The sharp telluric line features (e.g. located at 60 km s~^) are removed 
during the unsharp masking process. 

Speedy M ic to be one of the most active stars in the solar neigh- 
bourhood. |sjng^^3 found loK (^Lx/Lboi ) — 3 placing 
it in the saturated regime. More recently Mak^3 HqO^ reported 
log(Lx/Lboi) — —2 although this may be due to flaring as other- 
wise it would put Speedy Mic's X-ray luminosity an order of mag- 
nitude greater than any other star in its class. 

Prominences were first observed on Speedy Mic by I Jeffrie^ 

(l993). Two prominences were reported at radial distances of 
1.36 ± 0.09i?, and 2.96 ± 0.14i?.. Barnes et al. (2001) also re- 
ported four prominences at radial distances in the range 1.3 < i? < 
3.8J?.. 

In ij2|we present observations with uninterrupted and com- 
plete phase coverage of the Speedy Mic Hq line on consecutive 
nights. This strategy is combined with an improved technique for 
tracking prominences described in !j3| The relatively large number 
of resulting prominences allow us to examine the distribution of 
prominence heights and evolution timescales in We discuss the 
implications these results have on the nature of supporting mag- 
netic structures in S|5| Finally we present our conclusions in f|6| 



2 OBSERVATIONS 

Observations were made with the University College London 
Echelle Spectrograph (UCLES) at the Anglo Australian Telescope 
(AAT) on 2002 July 18, 19 &. 23. A slit width of 1.2 " was used, 
providing a mean resolution of ~ 45000 at Hq with coverage from 
4376A to 6892A in 47 orders on the 2Kx4K EEV CCD. 

Full phase coverage was obtained on both July 18 and 19 cor- 
responding to 100 and 97 spectra respectively. The typical signal- 
to-noise (S/N) achieved was 80 in a 300s exposure. On July 23 
time was lost due to bad weather and thus only resulted in 67 us- 
able spectra (~ 70% phase coverage) some of which have poorer 
S/N. 

Spectral extraction was perform ed using EC HOMOP, the 
echelle reduction package developed by Mills ( 1994) and its imple- 
menta tion of the optimal extraction algorithm developed bv lHomj 
il98d) . Further details of the extraction p rocess and a journal of ob- 
servations can be found in lBamedliool . Throughout this paper we 



M. Jardine 

Table 1. System parameters for Speedy Mic from lBamesI l2005h 



P [d] 0.38007 ± 0.00005 

[kms-i] -8.0 ±1.0 

^! sin j [kms^^] 132 ±2 

Axial inclination [deg] 70.0 ± 5 



shall be usi ng the system parameters found using the same dataset 
by Barnes ( 2005). These are summarised in TableQ We note that 
Iwolter et al. 1 , 2005i) find very similar parameters from data taken a 
few weeks later using the 8-m Very Large Telescope (VLT). 

We take the Ha order and normalise the continuum of each in- 
dividual spectrum. Any residual tilt in the spectrum is corrected by 
using a rigid spline fit after masking out the Ha line region. In Fig. 
Qan example observed Ha profile can be seen with a prominence 
absorption feature finishing its path across the stellar disc. The 
spectra are stacked together and phased according to their times 
of observation using midnight (UT) of the first night as the epoch 
of phase zero. 

The phased raw Ha time series are shown in Fig.|2| A visual 
examination of the time series reveals Speedy Mic to have a re- 
markably active prominence system. Many rapid transient absorp- 
tion features are seen to drift through the stellar Ha profile. 

The most noticeable feature however is a very dark region be- 
tween phases 0.8 < 4> < 1-0. Further examination suggests that 
the timeseries splits into two halves, with a lighter region between 
phases 0.2 < (f> < 0.7 and the rest of the timeseries being dom- 
inated by darker absorption features. We will briefly discuss these 
contrasting regions in The rest of this paper concerns the rapid 
transients seen in the Ha profile. 



3 PROMINENCE TRACKING 

The spatial distribution of prominences can be inferred from the 
timeseries. The phase at which any individual prominence transits 
the center of the stellar rotation profile can be easily identified. In 
addition to this by assuming that the prominences are in co-rotation 
with the surface of Speedy Mic (this assumption is validated by re- 
acquisition of the same prominence on subsequent nights) then we 
can calculate the axial distances (heights). This is due to the fact 
that the drift rates of the absorption transients through the stellar ro- 
tation profile are directly proportional to their radial diatnces from 
the rotation axis. Due to the relatively large heights of the promi- 
nences their orbital velocity sinusoids are well approximated by 
straight lines while they are transiting the stellar disc. 
The projected height is given by: 

■UJ V 

Rt Slusini 

where — RcosO is the distance of a prominence at latitude 9 
from the rotation axis, i) is the observed drift rate through the stellar 
rotation profile and is obtained directly from the phased timeseries. 
Q is the stellar angular velocity, i?* is the stellar radius and v sin i 
is the projected equatorial rotation speed. 

Previous authors have used a number of different methods 
to estimate the drift rates of the prominence absorption signa- 
tures through the line profile. When the S/N is poor, or there are 
gaps in the phase coverage of the prominence tracks, estimates 



Speedy Mic prominence system - 1 3 



2002 July 18 



2002 July 19 



2002 July 23 



T' 




-200 200 

Velocity (km/s) 



-200 200 

Velocity (km/s) 



'mm 



\mm 




I 
I 
I 

i I 
r I 







1 






■ 1 



-200 200 

Velocity (km/s) 



Figure 2. Raw Ha time series spectra for 2002 July 18, 19 & 23, with phase plotted against velocity. Dashed black lines show the v sini limits and the central 
dashed white line shows the radial velocity of Speedy Mic. The grey-scale runs from black at 0.80 times the local continuum level, to white at 1.02 times 
continuum. 



can be made directly by drawing a by-eye best fit line through 
the pro minence absorption signatures as was done in Barnes et al, 
i200lh . However the associated uncertainties were quoted as being 
~ 17?t. With more complete datasets the skew-mapping technique 
as devised in Collier Cameron & Robinson ( 1989a) and updated by 
IPonati & Collier Cameron 1 1997) can be used. In order to bring out 
just the ridge-lines of the absorption features an unsharp-mask is 
often performed on the raw Ha timeseries. Further details will be 
given in the next section. 

These techniques do not attempt to use a prior model for the 
prominences as they pass through the profile. The advantage of us- 
ing a model is that we can obtain a confidence assessment through 
a analysis. This describes how well the model fits the data and 
so is not as biased by strong areas of absorption but rather favours 
a consistent track across the profile. We describe the use of such 
a matched filter analysis technique with a Gaussian model in i|3.2l 
The first stage in the analysis however is to remove the atmospheric 
telluric lines from the timeseries spectra shown in Fig.|2| while at 
the same time bringing out the ridge-lines of the prominences. 



3.1 Unsharp masking 

Previous studies tracking prominences through the Ha profile have 
subtracted or divided by a mean profile in order to remove the tel- 
luric lines. As Fig.|2|shows however the Speedy Mic Ha profile is 
filled with many strong absorption signatures and it is difficult to 
establish a base reference level for the Ha line. The mean profile 
of the Ha line obtained by collapsing the timeseries in the temporal 
direction is therefore in strong absorption (see Fig.0. When the 
timeseries is divided by the mean profile we obtain spurious regions 



of relative emission, in which weak absorption transients can be 
seen. This is undesirable for further analysis as the absorption fea- 
tures should be below the continuum level. Another problem when 
using a mean profile is the assumption that the Ha profile is sym- 
metric, i.e. every transient absorption feature is symmetric around 
the centre of the stellar Ha profile. Again though just from look- 
ing at some of the absorption transients this not the case especially 
around the 0.8 < (j) < 1.0 region where an excess of absorption in 
the red is seen. 

In order to remove the telluric lines and to bring out the ridge 
lines of the absorption features we perform an unsharp mask in 
the temporal direction. We thus smooth the timeseries using a 
Gaussian-weighted running mean with a = 1500s. This corre- 
sponds to approximately half the duration of a typical absorption 
transient. We then divide the original timeseries by this smoothed 
version. The resulting timeseries spectra are shown in Fig.|4| This 
process has two main advantages. Firstly it removes the telluric 
lines very well as it does so on localised regions rather than the 
whole timeseries. This then allows for changes in the strengths 
of the telluric features as observations are made through varying 
amounts of the Earth's atmosphere during the night. Such a proce- 
dure also brings out the ridge lines of the transients making them 
more obvious. It is important to realise however that this is achieved 
because this technique essentially simplifies the timeseries, remov- 
ing the underlying slow variations of the Ha profile. Another point 
also worth remembering is that the absolute strengths of the absorp- 
tion (or indeed emission) with respect to the continuum level are 
not maintained. However this technique provides us with a higher 
contrast version of the timeseries which has brought out the rapid 
absorption transients that we are interested in. 
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3.2 Matched filter analysis tracking of prominence features 

In order to track the absorption features as they move across 
the V sin i profile of the star we use a matched filter analysis 
technique. This is based upon the spot traclcing technique of 
[collier Cameron et al. (2002) adapted for the prominence situa- 
tion. We model the absorption features as Gaussians of fixed 
width moving linearly with time through the Ha profile. We adopt 
a =25km s^^ which was found to be approximately the width of 
the weakest absorption transients. We then sample the timeseries 
at just over twice the number of spectra, a spacing of 0.004 in 
phase to ensure precise determination of prominence phases. At 
each phase the path through the profile of a potential prominence at 
that meridian crossing is considered for a range of possible promi- 
nence heights. Heights are considered from the stellar surface out 
to 7. 5i?* with a resolution of 0.0 15i?* to ensure adequate radial 
resolution for the low-lying prominences located well inside the 
co-rotation radius. Hence over a hundred thousand possible promi- 
nence tracks are tried throughout the timeseries. 

At each prominence height the model Gaussian is then opti- 
mally scaled to fit the data along the whole prominence track. The 
result is a single scale factor for the potential prominence moving 
through the v sin i profile at that particular phase and height, along 
with an associated x'^ for the fit. 

Initially all points that fell between the v sin i limits of the 
star were selected for Gaussian fitting, so that equal weight was 
given independent of position along the v sin i profile. This resem- 
bles a step function along the profile, with a value of one inside 
the V sin i range and zero outside. This approach was found to be 
an over-simplification when tested on actual absorption transients. 
The strength of a given absorption transient is determined by the 
physical properties of the cloud and its orientation with respect to 
our line of sight. If the cloud appears to transit the stellar disc at a 
high latitude then the extent of the absorption transient in the ro- 
tational profile will be reduced. The finite size of the cloud means 
that it takes time to move fully onto the disc. This combined with 
the velocity dispersion within the cloud itself and stellar limb dark- 
ening all change the strengths of the absorption features as they 
cross the disc. 

Furthermore for prominences at larger heights the path of the 
absorption signatures samples less spectra. Thus at the maximum 
height that we consider, the absorption signature samples only five 
spectra. This produces discontinuities in the resulting x^map, since 
as we increase the height we jump from considering n spectra to n-1 
spectra. 

In order to make a more realistic model for the absorption 
strength across the profile we parameterised all these factors into 
a single weighting function, shown in Fig.|3| This took the form of 
a Fermi-Dirac like function 



(l.0 + exp(a%^))' 

where the parameter a controls the nature of the fall-off away from 
the line centre and Aiiprom is an approximation to the internal 
prominence velocity dispersion. We approximated Awprom to be 
roughly twice the width of the fitted absorption features, a =25 
km s~^. The parameter a was changed to crudely optimise the fit 
to observed prominences passing through the profile, resulting in 
the profile shown in Fig.|3| This gives less weight to the extremes 
of the V sin i limits and it also produces smooth minima by re- 
moving the sharp and discontinuous jumps in the number of spectra 
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Figure 3. The weight given to the .scale factor of the Gaussian model as it 
crosses the stellar v sin i profile. Function is centi'ed on Speedy Mic's radial 
velocity, = 8km s"'^ 

sampled. Note that this function is of course only an approximation 
as the precise strength of any given feature is unique. 

Scale factors obtained through this procedure are either neg- 
ative or positive, corresponding to prominence absorption signa- 
tures or emission features respectively. It should be noted how- 
ever that most of the emission features are merely a product of the 
unsharp-mask process and are therefore not real. Just relying on 
the scale factor of the Gaussian is insufficient to reveal the promi- 
nence heights. This is because the scale factor alone will tend to 
overestimate the height of the prominences as it can join neigh- 
bouring prominences together. The heights we find are therefore 
based upon the cuts (at constant phase) through the resulting 
map at each phase considered. The resulting minimum in height 
is then compared to neighbouring phases (constant height) to estab- 
lish a local x^ minimum for a particular prominence. This results in 
much better determined minima and are used to produce the promi- 
nence tracks that are overlaid on the unsharp masked timeseries in 
Fig.E] Overlaying the tracks upon the timeseries in this manner al- 
lows quick visual confirmation that a suitable fit has been achieved. 

The uncertainty in the prominence heights can then be esti- 
mated from the curvature of the x^ slices. However due to the 
asymmetry in the x^ minima the reported uncertainties are mean 
estimates and almost certainty over-estimate the formal error on 
the prominence heights. The results of this improved tracking tech- 
nique on the Speedy Mic dataset are displayed in Table|2| 



4 PROPERTIES OF THE PROMINENCE SYSTEM 

The resulting prominence system is illustrated by the circular back 
projections displayed in Fig. |5| A total of ten prominences are 
found on the observable hemisphere of Speedy Mic on both nights 
with full phase coverage with a further five on the July 23 due to 
only partial phase coverage. For reference the most studied stellar 
prominence system, on AB Dor, is typically quoted as having 6 t o 8 
prominences at any one time (e.g. lCoUier Cameron et alll9 99h. It 
should be noted however that obtaining full and continuous phase 
coverage of a single rotation combined with applying the improved 
prominence tracking technique may affect this comparison. 

In total 25 prominences are found in the range 1.4 < Rc < 
3.6 stellar radii above the surface. This is the largest number of 
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Figure 4. As Fig.l2lbut after the unsharp masking process described in the text. The grey-scale now runs from black at 0.97 times the local continuum level to 
white at 1 .03 times continuum. Superimposed are the fits achieved by the prominence tracking technique. 



Table 2. The results of the prominence tracking analysis displayed with phases of meridian crossing and calculated heights for each day of observation. The 
quoted phases have a measuring uncertainty not greater than A<f> = 0.002 



Prominence 


Cycle 


July 18 
Phase -S- 


Cycle 


/((/)• 19 
Phase 


Cycle 


July 23 
Phase -S- 


A 


2 


0.102 


2.79 ±0.11 


4 


0.109 


3.18 ±0.16 


15 


0.116 


3.00 ±0.19 


B 


1 


0.219 


3.16 ±0.15 














C 








4 


0.267 


1.71 ±0.09 


14 


0.273 


3.09 ±0.17 


D 




0.343 


2.88 ±0.14 


4 


0.371 


3.31 ±0.26 








E 




0.430 


2.40 ±0.24 


4 


0.434 


3.00 ±0.42 


14 


0.407 


2.35 ±0.11 


F 




0.535 


3.28 ±0.18 


4 


0.547 


3.64 ±0.70 


14 


0.562 


3.39 ±0.38 


G 




0.609 


3.06 ± 0.34 


4 


0.620 


2.10 ±0.12 








H 




0.703 


2.82 ±0.37 


4 


0.726 


2.49 ±0.10 


14 


0.709 


3.24 ±0.17 


I 




0.813 


1.43 ±0.05 


3 


0.828 


2.87 ±0.40 








J 




0.871 


3.21 ±0.22 


3 


0.875 


3.48 ± 0.63 








K 








3 


0.919 


3.16 ±0.21 








L 


1 


0.969 


2.32 ±0.12 















prominences observed at any single epoch on any star thus far. 
Such a large number allows us to perform a statistical analysis of 
the prominence locations. The distribution of prominence heights 
is shown in the form of a histogram in Fig.|6| 



4.1 Prominence heights and the stellar co-rotation radius 

System parameters for Speedy Mic have been constantly i mproved 
over th e last decade. Using the same dataset as this paper iBamesI 
i2005l) were able to further constrain the period, stellar v sin i and 
axial inclination (see Table0. Such accurate parameters help us to 



identify Speedy Mic's co-rotation radius using: 

The co-rotation radius is essentially only as uncertain as the mass 
we adopt and even then is only proportional to the cube root of the 
stellar mass. This gives the absolute position of the co-rotation ra- 
dius however we need the co-rotation radius in terms of the Speedy 
Mic stellar radius for comparison with prominence height loca- 
tions. To calculate this we can use a simple dynamical argument 
assuming circular motion, using the observed v sin i, inclination 
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Figure 5. Circular back projections of prominence system. Pliase increases anti-clockwise, radial distances interior to the star have been suppressed and the 
dashed circle shows the location of the co-rotation radius. The grey scale represents the scale factor of the fitted Gaussian with black at -0.5 and white at -1-0.5. 
Black (negative) features show prominence absorption transients. Circular markers show the locations of the prominences from the x^analysis. The bottom 
right circle is a composite diagram of all prominences, those features identified on both July 18 and 19 have been joined by an aiTow and crosses are July 23 
data points. The dashed lines are radial and pass through the July 18 data point to aid the eye. 



and period. We perform a Monte Carlo error analysis using the val- 
ues and Gaussian-distributed la uncertainties from TableQ Note 
that for the uncertainty in the period we adopt the photometric pe- 
riod of 0.380 ± 0.004 days as given bv lCutispoto et alJil99"1) . As 
this uncertainty is based on a much longer observing run and so 
is probably more realistic. These parameters give a mean radius 
for Speedy Mic of i? = 1.06 ± 0.04Rq, w hich as one might ex- 
pect is considerably larger than the O.73i?0 jGravll992h Zero Age 
Main Sequence (ZAMS) radius that the K3V spectral type would 
suggest. Using Speedy Mic's determined radius and its associated 
error we can feed these back into another Monte Carlo trial to 
find the co-rotation radius in stellar radii. For the mass we adopt 
M« = 0.82 ± 0.08M(7). Th is covers a range of masses includ- 
ing the 0.74M,7. jGravlll992h K3V ZAMS and that considered by 
IWolteretalJ<200a) ^ using pre-main sequence evolutionary models 



of lSiessetai]i200d) of 0.9Mq. The location of the equatorial co- 
rotation radius is found to be Rk — 1.95 ± 0.07R, . 

The histogram of prominence heights shown in Fig.|6|clearly 
has a strong peak near 3_R«. This is considerably above the co- 
rotation radius, as illustrated in Fig. |6|by the hashed region. The 
actual mean of the prominence distribution is at 2.85 ± 0.54i?* 
with a median of 3.00R, which is less affected by the two promi- 
nences found below co-rotation. After the 3i?* peak the number 
of prominences falls off steeply with only two prominences found 
above 3.4J?». 



4.2 Evolution of prominence system 

This dataset provides us with the opportunity to study the evolution 
of the prominence system on two different timescales, that of con- 
secutive nights and then four nights later. In Table |2| we attempt to 
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Figure 6. Top panel: Distribution of 25 prominence lieights is sliown.Tlie 
solid line is tlie stellar surface and the hashed area shows the 1 a confidence 
interval on the location of the co-rotation radius. 



match prominences found on different nights based on the simple 
criterion of a phase difference less than = 0.03. This corre- 
sponds to approximately the extent in phase (the size) of the largest 
prominences. It is important, especially when considering consec- 
utive nights, to note what actual physical rotational cycle each part 
of the phased timeseries belongs to. The interval varies between 2 
and 3 rotations between the July 18 and July 19 and this is given in 
Table|2|for each identified feature. 

Using the above criterion to match the prominences we find 
that 8 of the 10 prominences observed on the July 18 have matching 
features on July 19. Furthermore that all 8 are observed at slightly 
later phases on the second night than on the first night. We also find 
that 6 of those 8 are observed to have moved outwards (increased in 
height) away from the star. These observations are shown in graph- 
ical form in the bottom-right panel of Fig.|5] 

It is worth briefly mentioning individual prominences in 
slightly more depth. Some prominences are easily identified on 
multiple nights. The strong absorption transient labelled as promi- 
nence A in Figs.|4|&|5|and Table|2|at phase <j) ~ 0.11 is present at 
all three epochs. Prominence F at </> ~ 0.54 is present at all three 
epochs also. In addition the faint absorption feature at ~ 0.87 is 
clearly visible on July 18 and 19. 

Prominence B at ~ 0.219 on July 18 seems to have disap- 
peared by the next night. Prominence G at iji ~ 0.61 is observed 
at a much lower height on July 19 than on July 18. We think this 
to be a newly formed prominence and in fact the feature observed 
on July 18 is still present but is much weaker than the new feature 
and is too blended in phase for the tracking technique to resolve the 
features unique minimum. 

Prominence H at </!> ~ 0.70 is faint on July 18 and at the start 
of observing on July 19 at phase ~ 0.73 a corresponding faint 
feature can just be seen finishing its track across the Ha profile. At 
the end of that night however a much stronger absorption feature is 
seen at a slightly later phase. We think this is most probably a new 
feature. This is interesting as it suggests a very rapid formation 
timescale. In less than 9 hours a new prominence has formed with 
a height ofR = 2.49i?.. 

Prominence I at </; ~ 0.82 is seen to experience a large in- 
crease in height between July 18 and July 19. This particular promi- 
nence signature is open to some interpretation. The fit on July 



18 is slightly ambiguous and the lower height is only marginally 
favoured over that of a prominence at a considerably greater height. 
Then the matching of prominence I with the feature on July 19 may 
just be co-coincidence. 

The observations on July 23 give us a valuable opportunity to 
study evolution on a longer timescale, with a gap greater than 10 
stellar rotations. Somewhat surprisingly, all five features found on 
July 23 have counterparts on July 19. The heights of the features 
appear fairly randomly distributed (and often in between) with re- 
spect to their July 18 and 19 counterparts. We therefore suggest that 
all five prominences have re-formed in the intervening time. 



5 DISCUSSION 

It is clear that the Speedy Mic prominence system is densely packed 
and relatively evenly spaced in longitude around the surface of the 
star. This high level of prominence activity is what might be ex- 
pected on this luminous X-ray ultra-rapid rotator. Interestingly we 
find that the prominences are concentrated at ~ 3-R*, well above 
the surface of the star and at twice the height of the co-rotation 
radius above the stella r surface. The traditional model of st ellar 
prominence systems by ICoUier Cameron & Robinsor] il989j) de- 
scribes prominences at a range of heights but with a concentration 
at the co-rotation radius. This dataset certainly challenges this as 
23 of the 25 prominences are found above the co-rotation radius 
with the peak of the prominence distribution at twice the height of 
the co-rotation radius above the stellar surface. This immediately 
requires large extended magnetic loops that go well beyond the co- 
rotation radius to support such structures. 

At first glance this would appear to be at odds with 
the high coronal densities measured with FU SE, Chandra 
and XMM-Newton of a range o f rapid rotators jDuBree.etalJ 
1 199.1"; 'Schriiver et al.' 1995; 'Brickhouse & Duoree' "199^; 
' Audard et al. 200 1 ; Mewe et al. 200 1 ; jyoung et al. 2001 ; 
Giidel et alj |2iroi|: |Sanz-Forcada. Brick house & Du pree 2003 ; 
iSanz-Forcada. Maggio & Micela .2003,) . These measurements 
suggest that the coronae of rapid rotators must be compact, since 
such high densities could not be confined at large distances above 
the stellar surface. The lack of significant rotation modulation 
in the X-ray emission can be explained if most of the emission 
comes from high-latitude regions that never pass behind the star 
as it rotates. This suggestion was given very striking support by 
the BeppoSax observations of two flares on AB Dor that showed 
no rotational modulation of the emission, despite the fact that 
the flares lasted for more than one rotation period of the star 
jMaggio et allbOOd) . Modelling of the decay phase of the flare 
suggested that the flaring regions were compact, and therefore 
must have been located close to the rotation pole. Emission from 
high latitudes on rapid rotators is not perhaps surprising given 
the very comm on appearance of high-latitude (even polar) spots 
on these stars I Stras smeied 1 1 996 ') . The relationship between the 
high-latitude flux and the coronal distribution of X-ray emission 
has been determined by extrapolating the surface magnetic field on 
AB Dor obtained from Zeeman-Doppler imaging and calculating 
the structure of the X-ray emission, assuming a simple hea t ing law 
ijardine. Collier Cameron & Donatil 120021: Ijardine et alj l2002h . 
The resulting X-ray emission is confined fairly close to the stellar 
surface and has a large component at high latitudes. By comparing 
the results of this technique with Chandra observations obtained 
simultan eously with the optic al spectra used for Zeeman-Doppler 
imaging jHussain et alj i2005h showed that not only the magnitude 



8 N.J. Dunstone, J.R. Barnes, A. Collier Cameron, M. Jardine 



and rotational modulation of the X-ray emission, but also the line 
shifts were consistent with the results of modelling based on field 
extrapolation. 

If the coronae of fairly rapid rotators such as AB Dor are com- 
pact, then it might be expected that at even higher rotation rates 
such as that for Speedy Mic the coronal extent might be even less. 
For the supersaturated star VXR45 Marino et al. (2003) found a 
large X-ray rotational modulation, showing that there is significant 
structure in the corona with clearly-defined bright and dark regions. 
This is what would be expected if the coro na had been centrifugally 
stripped as a result of the rapid rotation i Jardine & Unruhl|l999t 
[jardine 2004), leaving a corona that is very compact with only a 
patchy coverage of X-ray bright regions. In the case of Speedy Mic 
which appears to lie in the saturated part of the - period rela- 
tion it also appears unlikely that gas at coronal temperatures could 
be confined by the star's magnetic field at the heights at which 
the prominences are observed. The large number of prominences 
shows that the confining magnetic field must still possess a signifi- 
cant degree of complexity even out to these large heights. 

O ne recent theoretical model bv ljardine & Van Ballegooiienl 
<2005l) puts forward a new model for prominence formation that 
does not require a complex global magnetic field with large ex- 
tended loops in the X-ray corona supporting these prominences. 
Instead they propose that the prominences are supported by the 
re-connection of oppositely directed wind-bearing field lines over 
coronal helmet streamers. This would result in long thin highly 
curved loops extending out to large radii. These could then fill up 
with gas that was originally part of the wind a nd be dragged out- 
ward in radius. Ijardine & Van Ballegooiienl J2005i) find a maximum 
height for the summits of the coolest magnetic loops (i/m) in terms 
of the stellar co-rotation radius (Rk) given by: 




Using the Speedy Mic co-rotation radius Rk — 1.94/?* we 
find that ym ~ 3.35-R*. This would seem to be consistent with 
the o bserved prominence heights. Howev er it should be noted that 
while Ifedine & Van Ballegooiienl i2005l) illustrate that such mag- 
netic loops can achieve equilibrium they do not attempt to model 
the prominence structures themselves. 

The fall-off in the number of prominences detected beyond 
R ~ 3R, may suggest quite a different scenario to the one out- 
lined above. The traditional model of prominence support requires 
large magnetic loops in the X-ray corona. The strongest component 
though at large distances away from the star is the simple dipole, 
however, this could only provide effective support for prominences 
near the equatorial plane. Now if we assume that the observed 
prominences are actually in the equatorial plane then the highest 
prominence that would still transit the stellar disc assuming the in- 
clination of i = 70 ± 5° is at 7? ~ 3.73i?*. The highest detected 
prominence was feature F at 7? = 3.64i?, on July 19, so for Speedy 
Mic we cannot exclude the possibility that prominences are being 
supported in this manner. If this were the case, the highest promi- 
nences would only start to graze the stellar disc at high latitudes 
resulting in a smaller absorption signature. 

5.1 Evolution of the prominence system 

As illustrated in f|4|we have been able to identify prominences at 
similar longitudes on multiple epochs. On the timescale of con- 
secutive nights (2 to 3 rotations) we have witnessed the evolution, 



apparent formation and disappearance of several prominences. The 
fact that prominences are observed at similar phases more than 10 
rotations later is also interesting as it suggests considerable stability 
of the prominence supporting magnetic structures. 

Examination of the reported prominences heights and phases 
in Table|2|reveals a prominence system that appears to increase in 
height and lag behind in phase on the timescale of consecutive days. 
It should be noted though that the apparent changes in height of 
many of the prominences is still consistent, within the uncertainties, 
with a static prominence system. 

The observed phase lags are probably more accurate measure- 
ments. This is especially true of the narrow absorption features 
where the ridgeline of the absorption transients are distinct. The 
problem with the larger and often stronger absorption features is 
that they can be comprised of multiple components. This makes the 
confident re-identification of features on subsequent nights more 
difficult. 

The observed phase l ag between the epochs could be th e re- 
sult of a number of factors. IPonati & Collier CameronI ^99^ and 
lDonatietal]h999h attribute similar observations on AB Dor to the 
differential rotation of the prominence footprints. This then allows 
the latitude of the footprints to be located. In doing so though they 
assume that the supporting magnetic field lines are rigid and thus 
maintain prominence co-rotation. For Speedy Mic however we find 
that the phase lags between the consecutive days are generally of 
the order of A</> ~ 0.01. This results in prominence periods that 
are longer than the surface rotation rate at any latitude on the stel- 
lar surface and so cannot be explained by differential rotation. This 
may suggest that the stellar magnetic field is unable to enforce co- 
rotation. If we consider this in the context of prominences moving 
away from the star then the observed phase lag between consec- 
utive days could be the result of the prominences being dragged 
back in phase due to conservation of angular momentum. As such 
the phase lags observed between consecutive days have the poten- 
tial to inform us about the rigidity of the confining magnetic field. 

If the above model is correct then we need to explain why 
there is more noticeable phase lags between epochs on Speedy 
Mic than on AB Dor. This may simply be due to the location of 
the prominences relative to the co-rotation radius. We find that 
prominences group around a height which is twice that of the 
co-rotation radius above the stellar surface yet, as already dis- 
cussed, prominences on AB Dor are found at or just above the 
co-rotation radius. Therefore it maybe more difficult for the mag- 
netic field to enforce co-rotation of the prominence structures on 
Speedy Mic. Larger phase lags have however been observed be- 
fore on very high prominences on A B Dor. During observations in 
1994 lCollier Cameron et all <I999 h found that a prominence they 
observed moving radially outward from 3 or 4 stellar radii to 8 stel- 
lar radii lagged behind in phase. The authors attributed this to an 
attempt by the prominence to conserve angular momentum. 

Examining now the two prominences confidently identified 
on all three epochs, prominences A and F, we find that they both 
have relatively large phase lags between July 18 and 19 (two and 
three stellar rotations respectively). Then approximately ten rota- 
tions later we observe only a similar additional lag in phase on July 
23. This would suggest that the prominences have reformed at least 
once in the intervening ten rotations. This slower lag in phase may 
indeed be due to differential rotation of the prominence footprints. 
Barnes ( 2005) fitted a solar like differential rotation law to starspots 
observed on Speedy Mic of the form Q{9) — Qq — AQ,sm^{6), 
(where 9 is the stellar latitude and Slo the equatorial rotation veloc- 
ity). The values obtained were Q,o = 16.5361 ±0.0006 rad.day"^ 
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and AS7 = 0.033 ± 0.003 . We can use this to attempt to locate the 
latitudes of the prominence foot prints, as was done for AB Dor by 
iDonati & Collier Camer'or]ll99 7'). The observed phase lag over the 
ten rotations corresponds to a period for prominence A that would 
place it at a latitude of ~ 50° on the stellar surface. Similarly 
prominence F would correspond to a latitude of ~ 80°. The as- 
sociated uncertainties on these latitudes come from the uncertainty 
in the recurrence time and the uncertainty in the strength of the dif- 
ferential rotation. Resulting in an uncertainty of ~ 10°for promi- 
nence A and ~ 20° for prominence F. Thes e are si milar latitudes to 
those found bv lDonati & Collier Cameror] Jl997h . suggesting that 
the prominence footprints are located well above the stellar equator. 

Given the picture of prominence evolution that we have just 
described it seems that much care must be exercised in performing 
such an analysis of the prominence footprint locations. If the phase 
at which a prominence is observed is a function of its height then 
this must be taken into account. In order for a valid assessment of 
the phase lag due to differential rotation we need to observe the 
prominence at the same height on multiple epochs. It seems that 
this would require observing a prominence that has reformed at a 
similar phase and is at a similar stage in its height evolution. Due to 
our lack of understanding of the prominence formation mechanism 
we cannot be sure that we are really observing the same loop struc- 
tures. This is brought home by the data in Table |2| where we find 
that the two prominences (G and H) thought to have reformed be- 
tween July 18 and 19 were actually observed at slightly later phases 
than the original structures. This may point to a re-organisation of 
field at the prominence footprints such that the newly formed loops 
are not at exactly the same locations as the previous ones. If so 
then this would invalidate our analysis of the prominence footprint 
latitudes. 



6 CONCLUSIONS 

The improved prominence tracking technique allows for the fast, 
relatively automated detection and tracking of multiple promi- 
nences as they cross the stellar v sin i profile. This involved a 
matched filter analysis technique, fitting a gaussian model with 
a simple stellar limb-weighting function by a analysis. The 
Speedy Mic Ha timeseries has revealed one of the most (if not the 
most) densely packed and active prominence systems seen to date. 

We have shown that the distribution of prominence heights 
peaks considerably above the co-rotation radius. This provides fur- 
ther evidence against the traditional model of prominence support 
via large magnetic loops in the X-ray corona. Having discussed 
an alternative theoretical model of prominence support it would be 
very interesting to carry out Zeeman Doppler Imaging of Speedy 
Mic to map the surface magnetic field. In principle this could be 
done using SEMPOL at the AAT, but Speedy Mic is too faint to 
avoid excessive phase smearing. Unfortunately Speedy Mic has too 
southerly a declination for observation with the new echelle spec- 
tropolarimeter 'ESPaDOnS' (see Donati 2003) on the 4-m Canada- 
France-Hawaii Telescope. Ideally, an analagous high throughput 
spectropolarimeter at a 4-m or 8-m Southern hemisphere telescope 
is required. 

This study provides tentative evidence for a coherent picture 
of prominence evolution on Speedy Mic. Whether prominences are 
actually moving away from the star and lagging behind in phase 
as they attempt to conserve their angular momentum is still subject 
to interpretation. In order to more confidently distinguish between 
this effect and that of the differential rotation of the prominence 



footprints further observations will be necessary. These should fo- 
cus on tracking the evolution of individual prominence structures 
with the shortest time intervals possible. If we can obtain repeat 
phase coverage on three or more consecutive nights then we could 
hope to track multiple prominences through their entire evolution. 
Thus we could hope to define the relationship between prominence 
height and observed phase lag. This should give us an insight into 
the strength and nature of the magnetic support and further con- 
strain angular momentum loss through prominence ejection. 

While most of this paper has focussed on the rapid absorption 
transients interpreted as stellar prominences there are clearly other 
interesting things happening in the Ha line. The raw timeseries in 
Fig. m shows slower variations in the underlying Ha absorption. 
These seem to split the star into two longitude bands, for one half 
of the observed stellar rotation we see dark regions of absorption 
while for the half we see much less absorption. This is quite differ- 
ent from the Ha profile of AB Dor, where no such obvious contrast 
has been observed and as such requires further study. 
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